Magnetic and Thermodynamic Properties of the Three-Dimensional Periodic 

Anderson Hamiltonian 
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The three dimensional periodic Anderson model is studied with Quantum Monte Carlo. We find 
that the cross-over to the Kondo singlet regime is remarkably sharp at low temperatures, and that 
the behavior of magnetic correlations is consistently reflected in both the thermodynamics and the 
density of states. The abruptness of the transition suggests that energy changes associated with 
the screening of local moments by conduction electrons might be sufficient to drive large volume 
changes in systems where applied pressure tunes the ratio of interband hybridization to correlation 
energy. 

PACS numbers: 71.10.Fd, 71.10.Ht, 71.27.+a 
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The problem of localized, highly correlated electrons 
hybridizing with a conduction band is one of long- 
standing interest Q . Our understanding of the underly- 
ing physics has recently been increased through new an- 
alytic approaches |I],|2|], and numeric methods like Quan- 
tum Monte Carlo (QMC) These techniques have 
emphasized the connection between static magnetic prop- 
erties and the dynamic response like the density of states. 

However, what has been much less carefully explored 
by QMC is the link to thermodynamics. An intriguing 
problem for which a detailed understanding of the ther- 
modynamics is essential is the "volume-collapse" tran- 
sition in rare earth metals. This phenomenon occurs 
with the application of pressure to certain Lanthanidcs 
and gives rise to first order phase transitions with un- 
usually large volume changes (14% for Cerium and 9% 
for Praseodymium) |^|,||. Accompanying the change in 
volume is a change in the magnetism: On the expanded, 
highly correlated, side of the transition, the / electrons 
have well defined moments, while on the contracted, less 
correlated, side these moments disappear or are expected 
to disappear. The low- volume a phase of Ce is paramag- 
netic, as are the early actinides which are considered to 
be analogs for the collapsed rare earth phases . 

Even the qualitative origin of this phenomenon is still 
under debate. One suggestion is that the pressure- 
induced change in the ratio of the interaction strength 
to bandwidth gives rise to a Mott transition of the 4/ 
electrons accompanied by loss of magnetic order H . An 
alternate proposition is that the rapid change in the if— 
valence electron coupling leads to a "Kondo volume col- 
lapse" ||. In both cases, there are dramatic thermody- 
namic (e.g., pressure- volume) as well as magnetic signa- 
tures of the phenomenon. 

In this paper we will establish the connection between 
the thermodynamics and the magnetic properties of the 
symmetric periodic Anderson model (PAM) in three di- 
mensions. While previous efforts have focussed on the 



Anderson impurity model ||,[| , the capabilities of mod- 
ern massively parallel computers now make feasible rig- 
orous QMC calculations for the more realistic periodic 
model, which has so far received little attention in three 
dimensions. Our key results are: 

• The dependence of the singlet correlation function on 
the interband hybridization shows an increasingly sharp 
structure as the temperature is lowered, indicating a very 
rapid cross-over between a regime where the / sites have 
unscreened moments and one in which the moments are 
quenched by the conduction electrons. 

• A sharp thermodynamic feature exists at the same 
interband hybridization as this change in the singlet cor- 
relator. To analyze this, we introduce a new approach to 
the calculation of the free energy F, and show it obeys 
various analytic sum-rules. 

• The pressure difference at the transition inferred 
from F is reasonably consistent with experimental 
pressure-volume data on Ce, Pr, and Gd, given the ap- 
proximate representation of the electronic structure. 

The periodic Anderson Hamiltonian is 

ha k<7 

i 

+ ^ efUifg - fl y^X n ifv + n ida) ■ (1) 

We choose a simple cubic structure for which, 

£fc = — 2tdd [cos k x a + cos k y a + cos k z a] , 

Vk = —2tfd [cos k x a + cos k y a + cos k z a] , (2) 

where a is the lattice constant. The dispersion of Vk 
reflects our choice of near-neighbor (as opposed to on- 
site) hybridization of the / and d electrons |n| . Param- 
eter values and temperature T in this work are given in 
units of tdd- We take Uf — 6, consistent with the rare 
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earths for a reasonable choice of tdd = 1 eV, and explore 
a range of tfd and T values. Our choice of the "sym- 
metric" PAM dictates fx = tf = 0, and thus half-filling: 
(n>if) = (riid) = 1. QMC results for this model were ob- 
tained using the determinant algorithm H], which pro- 
vides an exact treatment (to within statistical errors and 
finite size effects) of the correlations. There is no "sign 
problem" for the symmetric PAM, allowing accurate sim- 
ulations at low temperatures. 

Figure 1 shows the temperature and t fd dependence of 
the singlet correlation function of near-neighbor sites i,j, 



C fd = (Sfi ■ S d j). 



(3) 



Here St 



and similarly for Sdy 



For weak interband hybridization, Cfd is small and the / 
moments are unscreened by the conduction electrons. At 
low temperature, a sharp change is seen to occur at tfd ~ 
0.6 to a phase where such screening is well established. 
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FIG. 1. The singlet correlation function (S f(i) ■ Sd(j)} as a 
function of f-d hybridization. As the temperature is lowered, 
there is an increasingly rapid switch from a small tfd regime 
where singlet correlations are absent to a large tfd regime 
where Kondo singlets are well formed. 

This sharp switch is also reflected in the energy E and 
free energy F. The difference AE(T) = E QMC (T) - 
-Eafhf(T) of the QMC calculations relative to antifer- 
romagnetic Hartree-Fock (AFHF) results at the same 
temperature (T = 0.08) is shown in Fig. 2. To get 
F = E - ST we fit pi the raw data for E QM c(T), 



Eqmc(T) = E + ^ c «e~ 



nA/T 



(4) 



The number of fitting parameters [Eq, c n , A) was taken 
to be about half of the number of data points. The en- 
tropy is then O] , 



S(T) = S + ^]>>(l + ^)e- 



^A/T 



(5) 



Fig. 3 shows a plot of the resulting free energy difference 
AF(T) = F QMC {T) - F AF kf{T). Independent fits {Eq, 
c„, A) were performed for each tfd, so that the smooth- 
ness of the resultant curves in Fig. 3 is one measure of the 
success of this procedure. Another is that our fit yields 
J2 n c n/nA to within ~ 3% of the expected value E| for 
tfd > 0.8. This sum is smaller by In 2 to within ~ 3% for 
tfd < 0.5, reflecting magnetic disorder of the spins below 
our lowest temperature (T = 0.08) in this regime, and 
consequent validity of the fit only for T > 0.08. 
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FIG. 2. The difference in energies between QMC and an- 
tiferromagnetic Hartree-Fock solutions. At small tfd, the 
AFHF energy accurately tracks the QMC. However, at in- 
termediate coupling the QMC results break away, reflecting 
the failure of HF to pick up the singlet correlations. A per- 
turbation approach (labelled E' 2 ') described in the text has 
some of the correct features seen in QMC. 

The crucial feature in Figs. 2 and 3 is the rapid change 
in slope at low temperatures of AE and AF near tfd — 
0.6. This behavior is hard to discern in the full thermo- 
dynamic functions whose variation with tfd is ~ 20 times 
larger than seen for these difference functions. It arises 
from the QMC results, and not from HF transitions, since 
the AFHF solution is stable throughout the tfd, T region 
plotted here. The size of the present slope change is not 
inconsistent with the volume collapse transitions, where 
one might view AF — min(Fi, F2) — F% with F\ and F% 
being free energy branches associated with the small t fd 
(large volume) and large tfd (small volume) phases, re- 
spectively. Given a volume dependence (f| of tf d ~ V~ 2 , 
the slope change is related to a pressure difference by 
VAP/2 = -{l/2)dAF/d\nV - dAF/d\nt Jd . Extrap- 
olations of experimental pressure-volume data ff^j into 
the two phase regions suggests VAP/2 ~ 0.4, 0.5, and 
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1.3 eV for Ce, Pr, and Gd, respectively. The low-T slope 
change in Figs. 2 and 3 is dAF/d\nt fd = 0.2-0.3 eV, 
which given the crudeness of the present representation 
of the rare earth valence electrons is reasonably consis- 
tent. 

The QMC calculations were carried out for a 4 3 - 
site lattice. As a systematic exploration of system size 
for these three-dimensional calculations would be pro- 
hibitive, we have used a second-order self-energy ap- 
proach to estimate the size effects, as well as to explore 
what analytic approximations might be more suitable 
than HF to capture the thermodynamics. The solid (4 3 ), 
dash-dot (6 3 ), and dash (8 3 ) curves in Fig. 2 labelled 
£( 2 ) were obtained from a finite-T version of the self- 
energy approach of Steiner et. al. jl5| . The Dyson equa- 
tion for the interacting Green's function matrix Gk is 
solved based on a second order (in Uf) expression for the 
self-energy E^ 2 '({G^}), determined from the paramag- 
netic HF result G k 0) . The trend in the AE^ curves for 
periodic clusters of 4 3 , 6 3 , and 8 3 sites suggests that finite 
size effects |l6| do not alter the qualitative physics and, 
indeed, move the position of the transition in the cor- 
rect direction for comparison with experiment, namely 
to higher values of dAF/dlntfd- 
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FIG. 3. The difference in free energies between QMC and 
antiferromagnetic Hartree-Fock solutions. At strong coupling 
(small tfd) the agreement in the free energy is good apart from 
an overall shift of Tin 2 associated with the tendency of HF to 
overestimate the magnetic order. As in Fig. 2, at intermediate 
coupling AF becomes sizeable. 

There is striking consistency between the singlet cor- 
relations in Fig. 1 and the energy and free energy differ- 
ences in Figs. 2,3. In all cases there is a rather abrupt 
switch in low-T behavior across tfd ~ 0.6, which anneals 
with increasing temperature. The anomalies are largely 
gone above T ~ 0.5, an upper bound for what might be 



a critical temperature in the present model. The actual 
critical temperature will reflect competition between ef- 
fects like these in AF and the volume dependence of a 
realistic generalization of Fafhf- An important term in 
AF is the QMC entropy, which reflects disordered spins 
for small tfd at the lowest temperature T — 0.08, in con- 
trast to both larger tfd values at this temperature, as 
well as the stable AFHF solution throughout the range 
plotted in Fig. 3, where the entropy is approximately 
minimal. Consequently, AF includes a —Tin 2 entropy 
term at small tfd, but not at large tfd, which serves to 
level out the AF curves as temperature is increased. 

Besides singlet formation, magnetic ordering of the lo- 
cal moments is a generic feature of the PAM. Indeed, our 
calculations of the /-/ structure factor suggest a strong 
tendency for the / moments to order antiferromagncti- 
cally at low temperatures with a maximal ordering tem- 
perature in the vicinity of tfd — 0.8. Further insight into 
the relation between AF, singlet formation, and the ther- 
modynamics can be obtained by computing the heat ca- 
pacity C(T) = dE(T)/dT §. We find a low-temperature 
peak similar to recent work on the two-dimensional Hub- 
bard model @, with an area J dTC(T)/T of In 2 at 
small tfd, which, however, washes out with decreasing 
area at large tfd- The peak has only minor impact on the 
slope change discussed above for AFqmc and AFqmc as 
functions of t fd C7j . 

A more complete picture of the PAM is given by the 
density of states, Nf(oj), which we obtain using the Max- 
imum Entropy method JL8| to perform the analytic con- 
tinuation of the imaginary time Greens function com- 
puted in QMC. The results for different tfd at fixed 
T = 0.2 are shown in Fig. 4. Nf(oj) evolves from a 
structure with upper and lower Hubbard bands separated 
by a gap Uf at small tfd to a regime where broadened 
remnants of these bands are still evident but additional 
resonant peaks characteristic of Kondo singlet formation 
have also developed. As tf d is increased, central reso- 
nances appear and are sharpest at tfd ~ 0.6, indicating 
the onset of singlet formation. Further increase in tfd 
enhances the weight in this central region at the expense 
of the Hubbard sidebands. 

The precise nature of the gap in the density of states at 
the Fermi surface, uj = 0, is still open to interpretation. 
For the half-filled, single band Hubbard Hamiltonian, 
N(u>) has a similar gap which evolves continuously from 
predominantly Mott-Hubbard character, for U >> W, 
to a Slater gap associated with antiferromagnetic order, 
for U << W . Similarly, the two-band model considered 
here has a Mott gap at small tfd, while the gap at larger 
tfd could originate either as a result of long range antifer- 
romagnetic order on the / sites, or, alternately, reflect a 
"coherence gap" associated with singlet formation. The 
competition between these two latter effects on N(uj) is 
well documented in a lower dimension [Q. Here, stud- 
ies of the /-/ correlation function show no signs of AF 
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long-range order at tfd = 0.6 and T = 0.2 > TncoI, which 
suggests these resonances signal singlet formation, not 
AFLRO. Analytic continuation of two particle Green's 
functions, like the magnetic susceptibility, will lend fur- 
ther insight into this question, but is very difficult and 
remains to be done. 
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FIG. 4. The /-band density of states for different f—d 
hybridization. For weak hybridization, there are peaks at 
iUf/2. These broaden with increasing tfd, and a Kondo res- 
onance develops. The curves for the two smallest tfd values 
have been reduced by a factor of two for display purposes. 

In this paper we have shown that there is a striking 
consistency between the location of sharp cross-overs in 
the singlet magnetic and thermodynamic properties of 
the three-dimensional periodic Anderson model. The / 
density of states shows a structure expected to arise from 
singlet correlations. Finally, estimates of the associated 
change in free energy are of the same order of magnitude 
as observed in the rare earth volume collapse transitions. 

Two important issues remain open. The first is the ex- 
tension to Hamiltonians with the full rare earth orbital 
complexity. Initial studies of how the Mott transition 
varies with band degeneracy in the Hubbard model, and 
other issues, already exist |l9|] within approximate nu- 
merical approaches like dynamical mean field theory §. 
The second, related, issue concerns band filling. Studies 
with many / orbitals will require working away from the 
symmetric point. 
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